ON THE Lp-SOLVABILITY OF HIGHER ORDER PARABOLIC 
AND ELLIPTIC SYSTEMS WITH BMO COEFFICIENTS 



HONGJIE DONG AND DOYOON KIM 



Abstract. We prove the solvability in Sobolev spaces for both divergence and 
non-divergence form higher order parabolic and elliptic systems in the whole 
space, on a half space, and on a bounded domain. The leading coefficients are 
assumed to be merely measurable in the time variable and have small mean 
oscillations with respect to the spatial variables in small balls or cylinders. 
For the proof, we develop a set of new techniques to produce mean oscillation 
estimates for systems on a half space. 
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1. Introduction 

The paper is devoted to the study of the Lp-theory of higher order parabolic and 
elliptic systems. More precisely, we expand the Lp-theory of higher order elliptic 
and parabolic systems to include a class of not necessarily continuous coefficients via 
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a unified approach for both, divergence type and non-divergence type systems in the 
whole space, on a half space, and on a bounded domain. The coefficients we consider 
are complex valued and, especially, the leading coefficients of parabolic systems are 
only measurable in the time variable and belong to the class of BMO (bounded 
mean oscillations) as functions of the spatial variables. The mean oscillations of 
the coefficients only need to be sufficiently small over small cylinders. 
To present the exact forms of systems, we let 

|Q|<m,|/3|<m |ct|^*Ti,|^|<m 

where m is a positive integer, 

and, for each a, (3, A°'^ = [A'^^ {t,x)]^ is an n x rt complex matrix- valued 
function. The involved functions are complex vector-valued functions, that is, 

The parabolic systems we study are 

ut + {-irLu = f, ut + {-ircu^ 

a|<m 

where the first one is in non-divergence form and the second one is in divergence 
form. The elliptic systems, non-divergence form and divergence form, respectively, 
are 

a|<m 

Whenever elliptic systems are considered, coefhcients, u, /, and are independent 
of t. When the domain is other than the whole space, we impose the homogeneous 
Dirichlet boundary condition. 

In the case of non-divergence type elliptic systems, we prove that, for a given 
/ G Lp(ri), there is a unique solution u G Wp"^{n) to the system Lu — f in H., where 
n is either the whole space R'^, the half space M^f. = {{xi, ■ ■ ■ ,Xd) G K"*,:!;! > 0}, 
or a bounded domain. We also prove the corresponding results for the other types 
of eUiptic and parabolic systems; see Section ||. 

As is well known, the key ingredient in establishing Lp-theory is apriori Lp- 
estimates of solutions to given systems. Largely, this is done in two steps. First, one 
establishes Lp-estimates for systems with 'simple' coefficients, for example, constant 
coefficients. Second, if the given system is in some sense close to systems with simple 
coefficients, one obtains the desired Lp-estimates by using a perturbation argument. 

The Lp-estimates for systems with constant coefficients, in many references, 
for example, ||], rely on the exact representation of solutions and the Calderon- 
Zygmund theorem. Another approach for such Lp-estimates is that of Campanato- 
Stampachia using Stampacchia's interpolation theorem (see |l9t). As to perturba- 
tion arguments, if the coefficients of given systems are uniformly continuous, the 
estimates are carried out by using the local closeness of the coefficients to constant 
coefficients in Loo norm. When the class of VMO (vanishing mean oscillations) co- 
efficients was first introduced, another perturbation argument was used in |0, ^, Q , 
where the continuity of coefficients is measured in the average sense, not in the 



HIGHER ORDER, SYSTEMS 



3 



pointwise sense, through a representation formula of solutions and the Coifman- 
Rochberg- Weiss commutator theorem. 

In this paper, in establishing the key Lp-estimates, we replace the first step, 
Lp-estimates of solutions to systems with simple coefficients, by mean oscillation 
estimates of solutions to the systems. Then for the second step we use a different 
perturbation argument, which is well suited to the mean oscillation estimates. For 
instance, if the system under consideration is elliptic in the form of Lu = / with 
constant coefficients in the whole space, then by the mean oscillation estimate of 
D^"^u we mean a pointwise estimate of the form 



/ \D^"''u- / D'^"'udy\dx 

J B^{xo) J BAxo) 



r(xo) J Br(xo) 

<Nk-^{1 |D2m^|2^^^ +Nj(-f \f\^dx^ (1.1) 



Co) / \J B^rixo) 



for all xo G K'^, r S (0,cx)), and n £ [ko,oo), where Br{xo) is a ball with center 
xq and radius r. Indeed, this implies the Lp-estimate of D^"^u by the well known 
Fefferman-Stein theorem on sharp functions and the Hardy-Littlewood maximal 
function theorem. But more importantly, this type of estimates well embraces 
the perturbation between the original systems and systems with simple coefficients 
when the coefficients have small mean oscillations over small balls or small parabolic 
cylinders. This approach was first introduced by Krylov to deal with second 

order elliptic and parabolic equations with VMO coefficients in the whole space, 
and is well explained in his book pfl] . 



Due to the well adaptiveness of estimates like (1.1) to the perturbation argument, 
our main effort in this paper focuses on obtaining mean oscillation estimates of 
systems with simple coefficients. Since in the parabolic case we allow coefficients 
to be only measurable in the time direction, the systems with simple coefficients in 
our case are naturally those with measurable coefficients depending only on t. 

For systems in the whole space, which corresponds to interior estimates, the mean 
oscillation estimates follow rather easily by adapting the techniques in ^ to 
higher order systems. However, differently from the arguments in p^ , we derive the 
non-divergence case as a corollary from the divergence case. Another noteworthy 
difference is that we prove the mean oscillation estimates not only for the highest 
order terms but also for the lowest order terms, so we are able to avoid the argument 
in deriving the Lp-estimates of solutions from those of the highest order terms, 
which is technically difficult in the case of higher order equations. 

For systems on a half space or on a bounded domain, which corresponds to 
boundary estimates, it is not possible to use the approach in ||2j, since the esti- 
mates developed there are only for equations in the whole space (interior estimates) . 
Thus here we develop a set of new techniques to produce mean oscillation estimates 
for systems on a half space. This is a new approach for boundary Lp-estimates, 
which is applicable to a wide class of equations or systems. To get these boundary 
mean oscillation estimates, as in the whole space case, we start with L2-estimates 
of systems on a half space. Although the L2-estimate for divergence type systems is 
well known under appropriate ellipticity or parabolicity conditions on the leading 



coefficients, our Theorem j.E regarding the L2-estimate for non-divergence type 



systems on a half space with coefficients measurable in time, as it alone, is a new 
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result to our best knowledge. In the proof we only use that of divergence type sys- 
tems and an interpolation argument. It is worth noting that L2-estimates for higher 
order elliptic equations and systems were obtained in [ |T8| , by using bootstrap 
arguments. For parabolic equations, however, in [ p^ the coefficients are assumed 
to be Holder continuous in the time variable since a semigroup method was used. 

From the L2-estimates, we derive the boundary mean oscillation estimates of 
some of highest order derivatives of solutions, precisely, -D™w in the case of diver- 
gence systems and D^u in the case of non-divergence systems, where x' denotes 
the last d—1 coordinates of a; = (xi, x') in W^. These estimates alone, however, are 
not sufficient for us to prove the main theorems. Because of this, we then consider 
a parabolic system with special coefficients, such that in a periodic pattern certain 
order normal derivatives of solutions to the system vanish on the boundary. This 
gives us the boundary mean oscillation estimates of D™u or D\"^u\ see Lemma 
3.4, Once we have all required mean oscillation estimates, we proceed as in [24 to 



the desired Lp-estimates using the perturbation argument, the details of which are 
illustrated for divergence type systems in the whole space; see Section ^. 

In the literature, for uniformly continuous coefficients, a rather complete Lp- 
theory can be found f or g eneral linear elliptic systems in ^] and for parabolic 
systems in [35, ^ If coefficients are in the class of VMO, non-divergence 



type higher order systems in the whole space have been investigated, for example, 
in |2^, ^ ^ , where leading coefficients of systems are either VMO with respect 
to all the variables or independent of the time variable. For divergence type higher 
order elliptic systems with VMO coefficients, we refer the reader to a recent in- 
teresting preprint [^^ in which the inhomogeneous Dirichlet problem on Lipschitz 
domains was studied. In all these papers, the method of singular integrals is used, 
so measurable coefficients are not allowed. 

Restricted to second order systems or equations, there are a relatively larger 
number of papers which can be compared to this paper. Non-divergence elliptic 
and parabolic equations on smooth domains with VMO coefficients were first stud- 
ied in H Q by using the technique of singular integrals. For further related 
results, we refer the reader to the book [^ and reference therein. The correspond- 
ing results for divergence elliptic equations were obtained in ^ by a similar 
technique. These results were later improved by the authors of [Ifin several papers 
for divergence type equations/systems without lower order terms on non-smooth 
domains by using a perturbation argument based on the maximal function theorem 
and a covering lemma (see Q for an extension to fourth order systems). An inter- 
esting question would be whether the methods in [^, ^ can be applied to equations 
with lower order terms or non-divergence form equations/systems. The methodol- 



ogy developed by Krylov in 24, g5| was later developed and extended in [|13[ for 
divergence and non-divergence systems in the whole space with the same class of 
coefficients, and in [|2[ |2^ for non-divergence parabolic and elliptic equations 
in the whole space with partially BMO coefficients for p > 2, and in for any 
p € (1, oo). In |ll|, this method was further adapted to divergence parabolic 
and elliptic equations/systems in the whole space with partially BMO coefficients. 
It is worth noting that in [0-jl4| and [^-[^ only interior mean oscillation esti- 
mates were derived. When dealing with equations and systems on a half space or 
on a bounded domain in |23, ESl O, H , the authors took full advantage of the facts 



that the coefficients are allowed to be merely measurable in one spatial direction 
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and the given systems are second order. Thus without using any boundary mean 
oscillation estimates developed here, the boundary Lp-estimates were derived from 
interior estimates as corohary type resuhs by using odd and even extension tech- 
niques. However, the extension techniques do not work for higher order equations 
or systems. This is the first paper in which the ideas in |2^, are adapted to 
boundary estimates, in both divergence and non-divergence cases. 

As noted above, the first critical step of the proof is the L2-estimates of systems 
with relatively simple coefficients under the ellipticity or parabolicity conditions 
on the leading coefficients. In this paper, we use so-called Legendre-Hadamard 
ellipticity condition, which is more general than the strong ellipticity condition 
considered, for example, in \29, |^, Nevertheless, it is still stronger than the 



uniform parabolicity condition in the sense of Petrovskii, which was used in jig , p3| , 
^ with more regularity assumptions on the leading coefficients. We shall discuss 
in details these conditions in Section ^ 

The organization of the paper is as follows. We introduce some notation and 
state the main results in the next section. The remaining part of the paper is 
divided into two parts. In the first part, we treat systems in the whole space. 
Section || and ^ are devoted to the L2-estimates and mean oscillation estimates 
for both divergence and non-divergence parabolic systems with simple coefficients. 
In Section || we complete the proofs of the Lp-solvability of systems in the whole 
space. The second part is the main part of the paper, in which we treat systems on 
a half space or on a bounded domain. In Section |^ we establish the L2-solvability of 
divergence and non-divergence parabolic systems with simple coefficients on a half 
space. Then in Section |7[ we obtain the boundary mean oscillation estimates of 
D^,u and D^Tu for divergence and non-divergence systems respectively. Section |^ 
is devoted to the estimates for a special type of systems. With these preparations, 
in Section |^ and |l^ we establish the Lp-solvability of both divergence and non- 
divergence parabolic systems on a half space and on a bounded domain. Finally, 
we discuss in Section some other ellipticity conditions used in the literature, and 
show how our results can be extended to systems under those conditions. 

2. Main results 

We first introduce some notation used throughout the paper. A point in R'^ is 
denoted by x = (xi,-- - ,Xd)- Whenever needed, we denote x by {xi,x') where 
x' G R''"^ A point in 



"1 = R X M'' = {{t, x) :teR,xe R'^} 

is denoted hy X — {t, x). For T £ (— oo, oo], set 

Ot = (-oo,T) X M'', £)+ = (-oo,T) X M^, 

where R"^ = {x = {xi, • • • , Xd) G R : xi > 0}. Especially, if T = oo, we have, for 
example, 0+ R x R'^. We also have 

Br{x) = {yeR'':\x-y\< r}, B^x') - {y' G R""^ : \x' - y'\ < r}, 

Qrit,x) = {t-r^''\t) X Br(x), QUt,x') = (t-r2'",f) X Blix'), 

Qt{t,x) = Qr{t,x)nOt,. 
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We denote 



Jn Jn 



For a function / on 2? C M + , we set 



= 1^ |" fit,x)dxdt ^ J f{t,x)dxdt, 



where I'D] is the d + l-dimensional Lebesgue measure of V. 

In order to state and prove our results on systems in Sobolev spaces, in addition 
to the well known spaces Lp and Wp, we introduce the following function spaces. 
As a solution space for non-divergence type parabolic equations, we use 

Wp^'2"((S',T) xn) = {u: ut,D"u £ Lp{{S,T) x f7),0 < |a| < 2m} 

equipped with its natural norm. Unless specified otherwise, in this paper D°'u(t, x) 
means the spatial derivative of u. For divergence type parabolic equations with 
= M'', we introduce 

VJ^{{S,T) X W^) = (1 - A)TWpi'2™((5'^T) X E-*) 
equipped with the norm 

l|u|k™((S,T)xR'i) = 11(1 - ^y~Aw^''^-^((S,T)x^'i)- 

Note that if we set 

Hp'"((5,r) X R'*) = (1 - A)^ip((S',T) X R"*), 

II-^I!h-'"((S,T)xR'') = 11(1 ^ ^) ~/IIlp((S,T)xR''), 

then 

ll"llw^((S,T)xR<') = ll"t|lH-'"((S,T)xR<i) + X! II -^""H ip((S,T) xR'') • 

|a| <m 

For a general 17, we set 

h;"((5, T)xn) = \f:f= D^u, u e Lpiis, t) X r!) 



\a.\<.m 



inf \ 



ll/llH-™((S,T)xn) - { Y ll/"llip((S,T)xn) '■ f - Y ^"-^^ \ ' 
l|Q|<7n |a|<m I 

and 

ni^iiS^T) xn)^{u:ute Hp"((S',T) x n),D"'u e Lp{{S,T) x 1]),0 < |a| < m}, 

l|w|k™((S,T)xO) = l|wt|lH-"((S,T)xn) + Y ll^""llip((S,T)xO)- 

I a I < m 

Let S, K > he two constants. Throughout the paper, we assume that all the 
coefficients are measurable, complex valued and bounded, 

'S'\ \a\ = \/3\=m, 
K, otherwise. 



I A"'' I < 



In addition, we impose the Legendre-Hadamard ellipticity on the leading coefficients 
(see, for instance, |l^). Here we call A°'^ the leading coefficients if |q;| — |/3| — 
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TO. All the other coefEcients are called lower-order coefficients. By the Legendre- 
Hadamard ellipticity we mean 

0''C£.'^A°''^{t,x)A >S\^\^"'\9\^ (2.1) 

\|a| = |/3|=m / 

for all {t,x) e M'^+i, ^ G R'', and 6 e C". Here we use 5R(/) to denote the real part 
off. 

Now we state our regularity assumption on the leading coefficients. Let 
osc^{A°''^,Qr{t,x)) ^ / \A"^{s,y)^ -f A"^ {s, z) dz\dy ds. 

Then we set 



^fl — ^^P ^^P 

(t,a;)GK''+i r<R. \a\ = \fi\=m 



We impose on the leading coefficients the small mean oscillation condition with 
a parameter p > 0, which will be specified later. 

Assumption 2.1 (p). There is a constant i?o G (0, 1] such that A"^^ < p. 

Contrary to non-divergence type systems where equations are defined almost 
everywhere, solutions to divergence type equations are understood in the weak 
sense. More precisely, for example, we say that u G T~{-p\^^{{S , T) x O), where 
I < p < oo, ft C R'', and -oo < S <T < oo, satisfies 

Ut + {-l)"'Cu + Xu= D"fa in {S,T) x n, 

\a\<m 

provided that 

/ f (^-ipfU+i-ir+\"\D°'ip-A°'^Df^ujdxds 

= (—1)'"' / / (f ■ fa dx ds + / u{S,x)(p{S,x) dx ~ / u{t,x)(p{t,x) dx 
Js Jn Jn Jn 



for every t G {S,T] and if = ((p\ • • • ,<p") G C°°((5,r) x Q) su ch that ipj t, ■) G 
C^{n) for aU tE[S,T]. If = -oo or T = oo, we take Lp G C^{{S,T) x fl) such 
that (p(— oo, •) = or ip{oo, •) = 0, respectively. 
We are now ready to present our main results. 

Theorem 2.2 (Divergence parabolic systems in the whole space). Let p G (1, oo), 
T G (—00,00] and fa G Lp{Ot) for \a\ < m. Then there exists a constant p = 
p{d,m,n,p,S) such that, under Assumption \2.i\ (p), the following hold true, 
(i) For any u G H^iOr) satisfying 

ut + {-l)"'Cu + Xu= Y ^"/" (2.2) 

|a| <m 

wc ho/VC 

Y X'~^\\D"uU^,^Or)<N Y >^^\\M\lAOt), 

I Q I < 7n I a I < m 

provided that A > Aq, where N and Ao > depend only on d, m, n, p, 6, K and 

Rq. 
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(ii) For any X > Ao, there exists a unique u e H™(Ot) satisfying (2.2). 

Theorem 2.3 (Non-divergence parabolic systems in the whole space). Let p £ 

(1, oo), T G (— oo, oo] and f S FMDt)- Then there exists a constant p = p(d, m,n,p, d) 
such that, under Assumption 2.1 (p), the following hold true. 

(i) For any u € Wp''^™{Ot) satisfying 

ut + i~iy''Lu + Xu = f inOr, (2.3) 

we have 

ll"t|lL,(0.)+ E >^''^\\D"uhAOT)<N\\f\\LAOr), 
\a\<2m 

provided that X > Xq, where N and Aq > depend only on d, m, n, p, 6, K and 

Rq. 

(ii) For any X > Xq, there exists a unique u G Wp''^"^{Ot) satisfying (2.3). 

Remark 2.4. We can also solve Cauchy problems for systems defined on (0, T) xW^ 
in divergence or non-divergence form. If the initial condition is zero, this is done 
by extending the original system to a system defined on (— oo,T) x K'' with the 
right-hand side being zero for t G (— oo,0). We deal with, in the same manner, 
Cauchy problems for the systems below defined on a half space or on a bounded 
domain. Note that in the case T < oo, by considering e^''^°~^^^*u instead of u we 
can take A = in the theorems above and below with the expense that N also 
depends on T. 

The next two theorems are about the boundary value problem of systems in 
divergence and non-divergence form on a half space — (— oo,T) x W^. 

Theorem 2.5 (Divergence parabolic systems on a half space). Let p G (l,oo), 
T G (—00,00] and fa G Lp{0^) for \a\ < m. Then there exists a constant p = 
p{d,m,n,p,d) such that, under Assumption \2.J\ (p), the following hold true. 

(i) For any u G W^{0^) satisfying 

ut + {~\TCu + Xu= D"fa inO+; 

\a\<m (2.4) 

It = Dim = ... = L>™"^u = ondpO^, 
where dpO^ ~ (— oo,T) x dWi^, we have 

E ^'"^II^""IImo+)<^ E (2-5) 

|Q|<m |a|<m 

provided that A > Aq, where N and Ao > depend only on d, m, n, p, 5, K and 

(ii) For any X > Xq, there exists a unique u G H™(Oj) satisfying (2.4). 

Theorem 2.6 (Non-divergence parabolic systems on a half space). Letp G (1, 00), 
T G (—00,00] and f G Lp{0 ^). Then there exists a constant p ~ p{d,m,n,p,S) 
such that, under Assumption \2.]\ (p), the following hold true, 
(i) For any u G Wp'^"^{0^) satisfying 

( Uf + (-l\"'Lu + Xu^ f in Ot; 

I , ^ ^ (2.6) 

1 w = L>iw = ... = L>5""^u = ondpO:^, 
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we have 

\a\<2m 

provided that A > Aq, where N and Aq > depend only on d, m, n, p, 5, K and 
Ro- 

(ii) For any A > Aq, there exists a unique u G Wp'^"^{0^) satisfying (2.6). 

Remark 2.7. By using a scaling argument, it is easy to see that we can choose Aq 
to be zero in the theorems above provided that £ or L has no lower-order terms 
and the leading coefficients depend only on t. 

Remark 2.8. In the above we presented the results only for parabolic systems. 
From those results we obtain easily the corresponding results for higher order elliptic 
systems in divergence form and non-divergence form. The key idea is viewing 
solutions to elliptic systems as steady state solutions to the corresponding parabolic 
systems. We refer the reader to Q and for details. To show the exact form 
of results for elliptic systems, we state below the cases for elliptic systems on a 
bounded domain. Theorem 2.12 and Theorem 2.13. 

Next we consider the solvability of systems in domains with the homogeneous 
Dirichlet boundary condition. For divergence systems, we assume the boundary dft 
of the domain n is locally the graph of a Lipschitz continuous function with a small 
Lipschitz constant. More precisely, we make the following assumption containing a 
parameter pi G (0, 1], which will be specified later. 

Assumption 2.9 (pi). There is a constant Ri £ (0, 1] such that, for any xq 6 Oft 

and r G (0, there exists a Lipschitz function (p: M.'^^^ — j> R such that 

Q n Br{xo) = {a; G Br{xo) : x'^ > 0(x')} 

and 

sup — < Pi 

in some coordinate system. 

Note that all domains satisfy this assumption for any pi > 0. Below we 
denote Qt (-oo, T) x Q, where fl C R'^. 

Theorem 2.10 (Divergence parabolic systems on a bounded domain). Let p G 

(1,00), T G (—00, 00] . Then there exist constants p = p{d, m, n,p, 5), pi = pi{d, m, n,p, (5, K, Rq) 
and Xq — Xo jd, m , n,p, 6, K, Rp, Ri) > 0, such that under Assumption |^. j| (p) and 
Assumption \2.!\ (pi) the following is true. For any fa G Lp(yi.T), \<y\ < Tn, and 
A > Aq, there is a unique solution u G ?^™(r2T) to 

+ (-l)™£u + Am = ^ D'^fa mVLr; 

|a| <m 

u = \Du\ = ... = ID^-^ul =0 on (-00, T) x dO., 
and we have 

J2 Ai-^p"u|U^(^,)<iv J2 a^II/.IImo.), 

|Q|<m |a|<m 

where N depends only on d, m, n, p, S, K, Rq and Ri. 
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Theorem 2.11 (Non-divergence parabolic systems on a bounded domain). Letp e 
(l,cx)), T G (— oOjCxd] and fl be a (7^™-!'! domain with the C'^"^^^'^ norm bounded 
by K . Then there exist constants p — p{d, m, 5) and Aq = \o{d, m, n,p, S, K, Rq) > 
0, such that under Assumption 2.1 (p) the following is true. For any f S Lp^flx) 
and A > Aq, there is a unique solution u G W^''^™' {Q,t) to 



ut + {-l)"'Lu + Xu = f in ^t; 
u = \Du\ = ... = \D"''-^u\ = on (-oo,r) x dVL, 



(2.7) 



have 



a I <2m 

where N depends only on d, m, n, p, S, K and Rq. 



As discussed in Remark 2.8, the theorems above have ehiptic analogies. We state 
the results below for elliptic systems on a bounded domain for future references. 

Theorem 2.12 (Divergence elliptic systems on a bounded domain). Letp G (1, oo). 
Then there exist constants p — p{d, m, n^p, S), pi = pi{d, m, n,p, S, K, Rq) and Xq = 
Xo{d,m^n,p, S, K, Ro, Ri) > 0, such that under Assumption 2.1 (p) and Assumption 
^.!\ (pi) the following is true. For any fa G Lp{Vl), \a\ < m and X > Xq, there is a 
unique solution u G W^{^) to 

\ol\ <m 

u = iZJul = ... = iD^-^ul = ondn, 
and we have 

a|<m |a|<m 

where N depends only on d, m, n, p, 5, K, Rq and Ri. 

Theorem 2.13 (Non-divergence elliptic systems on a bounded domain). Let p G 

(l,oo) and Q, be a (72m-i,i domain with the C'^"^^^'^ norm bounded by K. Then 
there exist constants p — p{d,m,n,p,d) and Aq — XQ{d,m,n,p, S, K, Rq) > 0, such 
that under Assumption 2.1 (p) the following is true. For any f G Lp[^) and X > Xq, 
there is a unique solution u G W^™(f2) to 

( Lu-f (-l)"Au = / infl; 

ju = = ... = iD^-^ui = on an, 

and we have 

\a\<2m 

where N depends only on d, m, n, p, 5, K and Rq. 
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Part 1. Systems in the whole space 

This part of the paper is devote d to the proo fs of the Lp-solvabihty of systems 
in the whole space, i.e., Theorem 2.2 and 2.3. In Section ^ we obtain several 
i2-estimates for systems with coefficients depending only on t. By using these 
estimates, in Section ^ we prove the mean oscillation estimates for systems with 



2.3 in 



the same class of coefficients. We complete the proofs of Theorem |2.2| and 
Section |[ 

3. L2-ESTIMATES FOR SYSTEMS WITH SIMPLE COEFFICIENTS IN THE WHOLE 

SPACE 

In this section we obtain L2-estimates of parabolic systems in divergence and 
non-divergence form when the coefficient matrices are measurable functio ns o f 
only the time variable satisfying the Legendre-Hadamard ellipticity condition (2.1). 
Even though our proofs are basic, we present them here for the sake of complete- 
ness. In particular, we derive the L2-estimate of systems in non-divergence form 
only using that of divergence type systems. Throughout the section we set 

Cou= J2 D°'{A°'^D^u), 

a I — I /3 1 — m 

where A"^ — A"^{t). Since A°'^ are independent of x G M'', we can write 

a I — I /3 1 — m 

Let C^{Ot) be the collection of infinitely differentiable functions defined on Ot 
vanishing for large 

Theorem 3.1. Let T £ (— cx),oo]. There exists N = N(d,n,m,d) such that, for 
any A > 0, 

J2 Al-i^P"zi|U,(0.) <7V J2 ^^WfJLAOr), (3.1) 
|a|<m |a|<m 

ifu e W^'iOr), U e L2{Ot), \a\ < m, and 

Ut + {-irCou + Xu^ D^U (3.2) 

in Ot- Furthermore, for A > and /„ £ L2{Ot), \ct\ < m, there exists a unique 
u e W^iOT) satisfying 

Proof. We assume A > 0. If A 0, the inequality (^^) holds trivially or we obtain 
Y <N Y \\fo.\\L,(Or) if /a = for \a\ < m 

\a\—m \a\—m 



using the inequality (3.1) for A > and letting A \ 0. 

Let us assume that the inequality (^.l[) is proved. Then due to the fact that 



ut = -(-1)'" D^iA^^D^i) -Xu+ Y D^fo.. 

\oi\ — \(3\—m |q I <7n 

we obtain \\u\\u^^Ot) < where Vxu ^ ut + {~1)"^Cqu + \u and 

N = N(d,n,ni,6,\). Then using the estimate, the method of continuity, and the 
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unique solvability of systems with coefficients A"^ — Saplnxn we prove the second 
assertion of the theorem. Therefore, we only need to prove the inequality (3.1). 
Moreover, since Vx is a bounded linear operator from W^{Ot) to ]HI^'"(C't), it 
suffices to concentrate on w G C^[Ot)- 

Multiply both sides of (3^) by u and integrate them on Ot- Then by integration 
by parts we have 

{u,Ut)o^ + {D''u,A'-PDPu)or+Hu,u)o^^ ^ (3.3) 

|a[ <m 

Note that 

JOt 

Here u is the Fourier transform of u in x. By the ellipticity condition we get 
sj \^\^"'\u\'^d^dt< [ ^(c^^u''A^\ d^dt. 



Also note that 

\u\'^(T,x)dx 



d_ 

dt 



\u\'^{t,x)dtdx = {u,Ut)oT + {ut,u)oj 



\u\^(T,x)dx > 0. 



'Ot 

Thus, if we denote the right-hand side of ( |3.3[ ) by /, we obtain 
5 



Since 



and 



/ \C\^"'\u\^dCdt + X{u.u)oT<'fiI < Y] \{D''u,U)ot\- 

JOt 

\{D^u,f^)oT\ < eX'^WD'^uWl^^a^^+Ne-'X-'^WUWl^^oT) 



for all £ > 0, the inequality (3.1) follows by using the interpolation inequalities and 
choosing an appropriate e. □ 

Theorem 3.2. Let T G (— oo, oo\. There exists N = N{d, n, m, 6) such that 
E ^'-^\\D"u\\L,(aT)<N\Wt + {-irCou + Xuh,^OT) 

\a\<2m 

for all X> and u £ W^'^"^{Ot)- Moreover, for X > and f G L2{Ot), there 
exists a unique u £ VF2^'^™(C't) satisfying 

Mt + (-l)'"£ou + Au-/ 

in Ot- 



Proof. As in the proof of Theorem 3.1, we only prove the estimate assuming that 
u G C^(p^). Let f ^ut + (-l)"/:ou + Au and write 

ut + {-l)"'D"{A°''^Df^u) + Xu^ f. (3.4) 



Then by Theorem 3.1 



A||u||l,(0.) <iV||/||L.(0.). 



(3.5) 
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Now by differentiating both sides of ( p.4[ ) m times with, respect to x we get 

{D"'u)t + (-l)™i:'"(A"^D^L>"u) + ALl^w = £>"/. 
This with Theorem |3.l| shows that 

^ \\D^D"^u\\L,(Or)<N\\fU,(o^y (3.6) 

\oi\—m 



Using ( p.q ), (3.6), and the interpolation inequalities, we obtain 

|Q|<2m 

Finally, observe that 

The theorem is proved. □ 

4. Mean oscillation estimates for systems in the whole space 
In this section we continue working on the operator 



u, 



\a\ — \p\—m 

where A'^^ = A">^{t). The main objective of this section is to obtain mean oscil- 
lation estimates for divergence type systems (Theorem ^!^) and for non-divergence 
type systems (Corollary 4.7) defined in the whole space. 

4.1. Some auxiliary results for systems in the whole space. First we prove 



the following localized version of Theorem 3.2 



Lemma 4.1. Let < r < i? < oo. Assume u G {Qr) and 

ut + (-irCou = f 

in Qr, where f G L2{Ot)- Then there exists a constant N = N{d,n,m,d) such 
that 

lktlU.(Q.) + P""«IIl.(Q.) < NWfhMu) + NiR - r)-^"luhMu)- (4-1) 
Furthermore, 

< iV|l/llL.(Q„) + N\\u\\L,iQ^), (4-2) 

where N = N(d, n,m,6, r, R) . 
Proof. Let 



1=1 



For each j = 0, 1, • • • , we take (j G C;f'(M'*+^) satisfying 
J 1 on Qr^ 

\0 on M'^+i\(-i?2™,i?,2™)xi3«^.^, 



C: 
and 



\D''Q\ < iV2'=^■(i^-r)-^ \iQ)t\ < N2^"'^{R-ry 
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where k ~ 0,1, - ■ ■ , 2m. Indeed, we can take (^j as follows. Let g{z) G C°°(M) be a 
function such that 

0<5<1, g(z) = 1 if z < 0, = if z > 1/2. 

Then set Q{t,x) = 4'jit)rij{x), where 

^,{t)^g{2^R^r)-\\t\^ ~R,)), 
r,,{x)=gi2^iR-r)-\\x\-R,)). 
Now we apply Theorem with A = to Qu G W2'^™{Oo), so that 

ll(0")tlU.(Oo) + P""(0")IIl.(o„) < N\\iQu)t + i-iy^Co{Qu)h,ioo) 

2m 

< mf\\L,(Qa)+miQ)tuh,^oo)+Nj2 p'o^""-'"IIl.(oo). (4.3) 

Using the properties of Q and interpolation inequalities (see, for instance, p6|), for 
each 1 < fc < 2m, we have 

<N2'^^R-r)-'^\\D'"^-''{C,+iu)h,^Oo) 
<£p2'"(0+iz.)|U,(Oo)+^22™^(i?-0"'"ll"llL.(Q„). (4.4) 
Furthermore, we have 

\\{Q)tuh,iOo) + WuD'^Qh^iOo) < A^22"^-(i?-r)-2"l«|U,(Q,). (4.5) 
Therefore, if we set 

= ll(0")t|U.(Oo) + P""(0")IU.(Oo), 

from (^!3|) , (0), and (|!|) we obtain 

/, < + iV||/|U,(Q,) + 7V22™J"(i? - r)-2™||y|U^(Q^). 

Multiply both sides by and make summations with respect to j to get 

CXD CXD OO OO 

Upon choosing, for example, e = 2^^™^^, the summations are finite, so from the 
above inequality we have 

p2"^(CoZi)||L.(Oo) + IIMtllL.(Oo) < ^II/IIl.(q«) +iV(i?-r)-2™l|u|U,(Q,). 
This proves the inequality ( |4.lD because the left-hand side of the above inequality 
is bigger than that of (^j|). Finally, the inequality (4.2) follows from (4.1) and the 
interpolation inequalities. □ 

In the sequel we denote u £ Wp°°{Qr{to,xo)), 1 < p < 00, (to, 2:0) € R'^+\ if 
D"u, D"ut € Lp{Qr{to,Xo)) for all multi-index a including a = (0, ■ • • ,0). 

Corollary 4.2. Let < r < R < (X) and u e W^^^iQa) satisfy 

Ut + i-iy^Cou^O (4.6) 
in Qr. Then for any multi-index 7, we have 

where N ~ N{d, n, m, 6, r, R, 7). 
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Proof. Note that 

in Qji. Hence it is enough to prove 



(4.7) 



Since u € W2 '' "^{Qr), this inequahty follows from ( [4.2[ ) if I7I < 2to, so assume that 
7I > 2m and 

D"'u = D^"'D'^u. 



Note that D^u is in wJ^'^^CQa) and satisfies (4.6). Thus applying (4.1) to the 
equation (4.6) with D^u in place of u we get 

where r < Rq < R- We repeat this process as many times as needed to get 

where |7o| < 2m and r < Ri < R. Then the inequality (4/7) for I7I > 2m follows 
from the same inequality for I7I < 2m (with i?i in place of r). □ 



(4.8) 



Lemma 4.3. //u £ W2'°°((54) satisfies ( [f.6[ ) m (54, i/ien 

sup|i:'u(i,x)| +sup|wt(t,x)| < N\\u\\l^(^q^), 
Qi Qi 

where N ^ N{d, n, m, 5) . 

Proof. Thanks to the fact that Ut = — (— l)'"£ow in Q4, it suffices to prove 

sup \D'^u{t,x)\<N\\u\\L^(Q^), 

for a multi-index 7. By the Sobolev embedding theorem 

\D^u{s,x)\^ds + N / \D''ut{s,x)\^ ds 

for each x G where D''u{t, x) is considered as a function of t e (—1, 0) for each 
fixed X Bi. On the other hand, again by the Sobolev embedding theorem there 
exists a positive number k such that 

sup \D'^u{s,x)\ < N\\D'^u{s,-)\\^,k(^g^^ 

for each s 6 (—1, 0), where D^u{s, x) is considered as a function of a: G _Bi for each 
fixed s G (—1, 0). We have the same inequality as above with D'^ut in place of D^u. 
Therefore, we obtain 



sup \D''uit,x)\'<Nj2 P''^^u|U,(Q,)+iV J2 P'^^^tlU.iQ, 



(t,2;)GQi 



i?|<fc 



|iJt<fc 



This together with Corollary 4.2 gives the inequality (4.i 

Lemma 4.4. Let A > and u e W2'°°(Q4) satisfy 

Mt + (-l)™>CoM + Au = 



□ 



(4.9) 
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in Qi. Then we have 

sup \D"'+^u{t, x)\ + sup \D"'ut{t, a;)| + sup \Du{t, x) 
Qi Qi Qi 



•A^sup|ut(i,x)| <7V^A^-^p''^u|U,(Q,), (4.10) 
Qi 



k=0 



where N = N{d, n, m, 5) . 



Proof. The case A = follows by Lemma 4.3 applied to D'^u since D'^u satisfies 
(4.(;). For the case A > 0, we follow an idea by S. Agmon. Consider 

C(?/) = cos(A^y) + sin{\^y). 

Note that 

{-irDi^ay)^\ay), m = i, \D^m\ = ^^- 

Denote by (t, z) = {t, x, y) a point in where z = (x, y) e K''+-'^, and set 

u(t, z) = u(t, x)C(y), Qr = (-r^'", 0) X {|z| < r, z e M'^+i}. 
Since u satisfies ([4.9|) , ii satisfy 



in Q4. Upon applying the inequality (4.10) with A = just proved above, we get 
sup \DZ'+^u{t, z)\ + sup \D'^Ut{i. z)\ + sup \D'^D^u{t, z)\ 



Qi 



Since, for example. 



+ sup \D:;Mt, ^)l < A^P"^i|L,(Q,). (4.11) 



sup A2 |D^u(t,x)| < sup \D'^D^u{t,z)\, 

it,x)£Qi {t,z)eQi 

the left-hand side of (4.11) is bigger than that of ( 4.10| ). On the other hand, D"^u 
is a linear combination of terms like 

Ik It, Ik 1 ] 

\2-2^ cos{X^y)D^u{t,x), X2^27^ sm{\^y)D^u{t,x), fc = 0, • • • , m. 

Thus we see that the right-hand side of ( [4.1lD is bounded by that of ( 4.1C ). The 
lemma is proved. □ 



Recall that we denote by X a point in 



prf+i 



Lemma 4.5. Letr G (0, 00 ), k e [4, ex)), A > 0, and Xq = {tQ,xo) G M'^^"'^. Assume 
u £ ^™;oc(IR'^~'"^) satisfies (O) in QK.r{Xo). Then for any a, \a\ — m, we have 

{\D"U - (i?"«)Q.(Xo)|)Q^(^„) + (1^. - (")q.(Xo)I)q„(;,„) 

771 

<N^-'Y.xi'^{\D'uf)l^^^^^^, (4.12) 

k=0 

where N = N{d, n, m, S) > 0. 
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Proof. Let us prove the inequality ( 1.12 ) when Xq — (0, 0). This with a translation 

of the coordinates proves the inequality for general Xq e R'^+^. 

Since the standard mollification of u with respect to x satisfies (4.9) in a little 
bit smaller cylinder than Q^r, we assume that Z3"m G L2{QKr) for all multi-index 
a. Furthermore, (4.9) implies that D'^ut G L2{QK.r) if D"u G L2{QK,r) for all a. 



rl,00 



Therefore, without loss of generality we assume that u G 

Due to a scaling argument (for instance, see the proof of Lemma |7.5| ), it suffices 
to deal with the case r = 4/k. Observe that, for example, 

(|£i"u - (D"u)qJ)o < NrsvLp\D"+'^u{t,x)\+ Nr sup \D''ut{t,x)\. 

Qi Qi 



By L emma [4.4| , the right-hand side of the above inequality is bounded by that of 
( 4.12| ) (recall r — 4k~^). The lemma is proved. □ 

4.2. Mean oscillation estimates for systems in the whole space. In the 

next theorem, we prove a mean oscillation estimate for divergence form systems 
with simple coefficients in the whole space. 

Theorem 4.6. Let r G (0,cxd), k g [8,oo), A > 0, = (to, So) e W^+^, and 
la G L2,ioc(K''+^); \a.\ < m. Assume that u G H^iocO^'^^^) satisfies 



ut + i-l)"" Cqu + Xu 



\a\ <m 



in QKriXo). Then for any a, \a\ — m, we have 

{\D'^U - (i?"«)Q^(X„)|)Q^(^„) + {\U - («)q.(Xo)|)q^(;,„) 



<7Vk-i^A5- 

fc=0 



Nn-^+i A^ 

I Q I < m 



(4.13) 



where N = N{d^ n, m, 5) > 0. 



Proof. We take, for the sake of simplicity, Xq = (0,0). As mentioned earlier, a 
translation gives the result for general Xq. 

Take an infinitely differentiable function C, defined on M'^+^ such that 

C = l on Q,,/2, C = outside (-(Atr)^'", («r)2™) x 



By Theorem 3.1, for A > 0, there exists a unique solution w G ^^{Ooo) to the 
equation 

Wt + {~l)^CQW + \w= D^iCM (4.14) 



|a| <m 



1. Let w u - w. Then the function v G H^iocO^"^^^) satisfies 
vt + (-l)™>Cot^ + Aw = in Q^rl2- 



By Lemma 4.5 (note that k/2 > 4) applied to w, we have 

m 



fc=0 



/2 

(4.15) 
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Next we estimate w, which is the unique solution to the equation (4.14) consid- 



ered on Oo- By Theorem 3.1, we have 



i2(Oo)- 



: I <m 



I Q I < m 



In particular, 



\a \ <m 



m ^ 



(4.16) 
(4.17) 



k=0 



\a\<m 



Now we are ready to prove (4.13). Since 

i\D'^u^iD"u)Q,.\)Q^<2i\D'^u~c\)Q^ 

for any constant c, by taking c — {D°'v)q^ and repeating the same argument for 
the second term, we bound the left-hand side of (4.13) by a constant times 

{\D'^u-iD'^v)Q,.\)Q^+xHu~{v)Q^)Q^, 

which is, due to the fact that u ~ w + v, controlled by 

Using ( [1.15| ) and ( fl.l6| ), we see that the above is less than 



k=0 



a I < m 



Finally, we use the fact that v — u ~ w and ( 4.17 ) to prove that the terms above 
are not greater than the right-hand side of ( 4.13 ). □ 

Next we consider the corresponding mean oscillation estimate for non-divergence 
type systems in the whole space. 

Corollary 4.7. Let r e (0,cx)), k € [8,oo), X > 0, Xo € R x R'^ , and f e 

^2,ioc(K'^+^), |a| < m. Assume that u e M^2jo™(K'*+^) satisfies 

ut + {~lY'Cou + Xu = f 
in QKriXo). Then for any a, \a\ = 2m, we have 

(|Z?"U - (^"")q.(Xo)|)q^(Xo) + ^ (I" - (")q.(^o)I)q,(Xo) 



k=0 



where N — N{d^ n, m, 6) > 0. 

Proof. Again let Xq = (0,0) for simplicity. By Theorem 4.6, it follows that (after 
multiplying both sides by A^) 

m 

A(k-(^^)Qj)Q. <iV^-'E^'"^(l^'"l')L+^'^"^^(l/l')L- (4-18) 

k=0 
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Differentiate m times both sides of tlie system witli respect to x to get 



By Tlieorem 4.6 applied to D'^u in place of u, 



fe=0 



where I7I — m. This combined with (4.18) gives the inequality in the corollary. □ 



5. Lp-ESTIMATES FOR SYSTEMS IN THE WHOLE SPACE 

In this section, we use the mean oscillation estimates obtained in the previous 
section to prove Theorem 2.2 and 2.3. 

Let Q= {Qrit.x) : {t,x) eM''+\r e (0,cx))}. For a function 5 defined on M'^+i, 
we denote its (parabolic) maximal and sharp function, respectively, by 



Mg{t, x) 
9*{t,x) 



sup + \g{s,y)\dyds, 
QeQ:(t,x)eQ J Q 

sup / \g{s,y) - {g)Q\dyds. 
QeQ:{t,x)eQ J Q 



Then 



\\g\\L,<N\\g*\\L,, IIMglU, <iV|l.g||L,, 



a g € Lp, where 1 < p < 00 and N — N{d,p). As is well known, the first inequality 
above is due to the Fefferman- Stein theorem on sharp functions and the second one 
is the Hardy-Littlewood maximal function theorem. 
We use the idea of freezing the coefficients to obtain 



Lemma 5.1. Let C be the operator in Theorem \2j_^. Suppose the lower-order co- 
efficients of C are all zero. Let fi,v € (1, 00), -\- l/v = 1, and A, i? G (0, 00). 
Assume u G {M.'^'^^) vanishing outside Qb, and 

ut + [-ircu+\u= J2 ^"/"' 

I Q I < m 

where fa G L2joc(^'^^^) ■ Then there exists a constant N = N{d,m,n,S, ij) such 
that for any a, \a\ — m, r E (0, 00), k > 8, and G M'^^"'^, we have 

d^""- (^"")q.(Xo)I)q.(Xo) + A^ (I" - («)q.(^o)I)q,(X„) 
m 

k=0 



. |a|<m 



, (A#)A(|£)™u|2^)^^ 
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Proof. Let K > 8 and r G (0, oo). Fix a ?/ e M'* and set CyU = A°'^{t, y)D'^D^u{t, x). 
Then we have 

I a I < m 

where 

U = U + {-ir E iA"^it,y)~A'^^{t,x))DPu. 

|/3|=m 

It foUows from Theorem 4.6 that 

{\D"u - (i?"«)Q.(x„)|)Q„(^„) + {\u - {u)QM\)Q^iXo) 

<^--'E^^"^(i^'"i')L(^o) + ^-"^^ E ^^'^(i/"i')L(xo)- (5-1) 

k—0 \a\<m 



Note that 



/ \fa?dxdt<N I \fafdxdt + Nly, (5.2) 



where, for |a| — to, 



ly^ liA^^it^y) - A"l^{t,x))Df'ufdxdt. 

Denote B to be Bf^rixo) if Kr < R, or to be Eft otherwise; denote Q to be QK,r{tQ, xq) 
if Kr < R, or to be Qr otherwise. Now we take average of ly with respect to y in 
B. Since u = outside Qn, by the Holder incquahty we get 

Iydy= / / \iA''^{t,y) - A''^{t,x))D^uf dxdtdy 

< i ( f \A"P{t,y)-A"P{t,x)\^A dy( [ \D"^u\A ' , 

J B \Jq / \JQnriXo)nQR ) 

where, by the bomrdedness of as well as the definitions of osc^: and A^, the 
integral over B in the last term above is bounded by a constant times 

|^"^(t,y)-v4"^(i,x)|V'dy< f / / \A^^{t,y)~A'^f'[t,x)\dxdtdX 

B \JQ J \J B JQ J 

< N (|g|osc,(A"^, Q)) ^ < iV ) " . 

This together with (|5.l|) and (|5.2|) completes the proof of the lemma. □ 



Proof of Theorem 2. St. Due to the method of continuity, it suffices to obtain an 
apriori estimate. By moving all the lower-order terms to the right-hand side and 
taking a sufficient large A, we may assume that all the lower-order coefficients are 
zero. 

Case 1: p € (2, oo). First we suppose T — oo and u G C^[QjiQ). Choose a /i > 1 
such that 2^ < p. Under these assumptions, from Lemma |5.l| we easily deduce 

k=0 
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+Nk 




for any a, |a| = m, r ^ (0,oo), k > 8, and Xq G M''+^. This together with the 
interpolation inequality, the Fefferman-Stein theorem and the Hardy-Littlewood 
maximal function theorem yields 



fe=0 



< Nin-i + pA) ^ A^-^^ \\D''u\\ 



Nk 



k=0 



E 

I q| <r) 



A^-^||/a||. 



(5.3) 



Now we can choose k sufficiently large and p sufficiently small in (5.3) to get the 
desired estimate. A standard partition of the unity enables us to remove the re- 
striction that u G C^{QRg). The extension to the case T € (— c)o,+oo] is by now 
standard; see, for instance, 24 . We omit the details. 



Case 2: p £ (1,2). Since the system is in divergence form, this case follows 
immediately from the previous case by using the duality argument. 

Finally the case p = 2 is obtained by an interpolation argument. □ 



In a similar way, from Corollary 4.7 we get the following counterpart of Lemma 
5.1 for non-divergence systems. 



Lemma 5.2. Let L be the operator in Theorem 2^. Suppose the lower-order co- 
ejficients of L are all zero. Let p,v € (l,c»), l/p + l/v = I, and A,i? G (0,oo). 
Assume u G (R'^'^^) vanishing outside Qr and 

ut^ {-!)"' Lu + f, 

where f G L2^ioc{^'^^^) ■ Then there exists a constant N — N{d, m, n, S, p) such that 
for any a, \a\ = 2m, r G (0, oo), k > 8, and Xq G M'^^"'^, we have 

- (^"")q.(Xo)|)q^(Xo) + a (I" - (")q.(^o)I)q,(Xo) 



fe=o 



Proof of Theorem 2.S. As in the proof of Theorem 2.2, it suffices to prove the apri 



ori estimate for T — oo. 

Case 1: p G (2,oo). We only need to consider the case when u G C^{QRg), 
since the general case follows from a partition of the unity. The proof of this case 
is almost the same as that of Theorem 2.2, by using Lemma instead of Lemma 
5.1. So we omit it. 



Case 2: p G (1,2]. Note that here we cannot use the duality argument directly. 
From Case 1 and Remark 2_/7, we already have the W^''^™ solvability of 

in the whole space for any q G (2, oo) and A > 0. For this system, since A"^ are 
measurable function of time only we can make use of the duality argument, which 
yields the solvability of the same equation for any q G (1,2). Fix aq = (l+p)/2. 
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Now we can repeat the arguments in the previous section with q in place of 2, and 



get the estimate in Lemma 5.2 with q in place of 2. Finally, following the proof of 



Case 1 completes the proof of this case. □ 

Part 2. Systems on a half space or a bounded domain 

This is the most novel part of the paper. The objective of this part is to establish 
the Lp-solvability of parabolic systems on a half space or on a domain. 

In the next section, we prove the L2-estimates for systems with coefficients mea- 
surable in t on a half space. Relying on these L2-estimates, in Section |7| we are 
able to derive mean oscillation estimates of some partial derivatives of solutions to 
systems on a half space. These estimates alone are not sufhcient for our purpose. 
So in Section |^ we consider a certain system with special coefficients. Combining 
the results in Sectio n fi\ a nd | together enables us to prove the Lp-solvability on a 
half space (Theorem 2^ , ^ ) . Section |l^ is devoted to the proofs of the bounded 



domain cases (Theorem 2.1C , 2.11 ). Finally we give several remarks about other 
ellipticity conditions. 

6. L2-ESTIMATES FOR SYSTEMS WITH SIMPLE COEFFICIENTS ON A HALF SPACE 

In this section, we prove the L2-estimate for systems on a half space. We again 
consider 

|q 1^1/3 |—m |q 1^1/3 |—m 

where A"^ = A°''^{t). RecaU that = (-00, T) x R^. In the divergence case 
(Theorem |6.1[), the proof is rather standard. However, in the case of non-divergence 



systems (Theorem S.6), the proof is much more involved. To the best of our knowl- 
edge. Theorem 5.6 



is new for higher order parabolic equations and systems with 
measurable coefficients depending only on t. 

6.1. Divergence case. Throughout the paper, we use the notation D™;u to in- 
dicate one of D°'u, where a — (ai,-- - ,0;^), ai = 0, and \a\ — m. Sometimes, 
depending on the context, D™;u means the whole collection of D^u, \a\ = to, 
ai — 0. Similar to C^{Ot), we denote by C^(O^) the collection of infinitely 
differentiable functions defined on vanishing for large |(i,a;)| G O^. 

Theorem 6.1. Let T G (—00,00] and fa G L2{0^). There exists a constant 
N — N(d, n, TO, S) such that 



1^ V-^P"^IL,(o+)<^ 1^ (6.1) 

|a|<m |a|<m 

for any A > and u G TL^{0^) satisfying 

u(t, 0, x') = • • ■ = D-^^^u{t, 0, x') = (6.2) 
on (-00, T) X W^-'^ and 

ut + {-irC^u + Xu^ Y ^"/a (6.3) 

I a I < m 

in . Furthermore, for A > and G L2{0^), \a\ < m, there exists a unique 



u G W^{0+) satisfying (|6j) in 0+ and (|6^ on (-00, T) x 
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Proof. As in the proof of Theorem 3.1 , we consider only the case A > 0. We 
foUow the hnes of the proof of Theorem 3.1. One noteworthy fact is that, because 
u e W^{0^) satisfies we have 



where the function u on the right-hand side is viewed as an extension of u to Ot 
so that it is zero on Ot \ . Similarly, 



JOt 



where u is the Fourier transform of the extension. 



□ 



Remark 6.2. Theorem 



6.1 



can be extended to systems in a cylindrical domain 
fix, where f2 is a bounded Lipschitz domain. For small A > 0, we have a better 
estimate than (3.1). Indeed, from the proof above, we get 



|q I <m 



By using the Poincare inequality, 

IWh^inr) < N\\Du\\L,in^) < N\\DML,inr) < - < iVp™^.|U,(o,). 
Thus, wc conclude 



I a I < m 



fe=0 



Note that in this case, the solvability also holds for A = 0. 

6.2. Non- divergence case. Let us introduce some additional notation. Let t G N 
and {ci, ■ • • , C2r} be the solution to the following system: 



fe=i 



2t- 1. 



(6.4) 



For a function w defined on R5j_, set 

w{xi, x') if xi > 

SrW = < 



Cfew(— — xi, a:') otherwise 
. fe=i 



Note that SrW e C'^^-^{W) ii w e C°°(]R« ). Indeed, by (|6J) 

Di (^^Ckw{-^xi,x')^ CkDiw{0,x') = DiwiO,x') 

for j = 0,-- - ,2t-1. 

We remark that similar extension operators were used in [^Sj and in the 
study of elliptic systems. We will use the following interpolation estimate. 
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Proposition 6.3. Let 1 < p < oo and u e W^{W\^). For any e > 0, there exists 
N — N(d, n, m, p, e) such that 

d 

where /c = 0, 1, • • • , m — 1. 

Proof. Without loss of generality we assume that u G C,5"(I'^+)- Let u — ErU. For 
a suffi-ciently large r, the extension u is in Wp^{W^) and satisfies X]^=i Dj"^w = 
Eti^r/^'^here 



U, 
2t 



Here Ck are appropriate constants. Observe that 

d 

where the second inequality is due to the Lp-estimate of elliptic systems in the 
whole space (see Remark p^ ) and N = N{d,n,m,p). By replacing u{xi,x') by 
u{eiXi, x') in the above inequality we have 

d 



2m „ 



J=2 



The proposition is proved. 

Lemma 6.4. Lei T G (— oo,cx)]. There exists N = N{d,n,m,S) such that 

\a\—m 

for all X > and u G W2''^™{0^) satisfying 

u(t, 0, x') = • ■ ■ = D'{'-'^u{t, 0, x') = 0. 

on (-oo,r) X R*^-!. 
Proof. Define 

in O^. Then by Theorem 



□ 



(6.5) 



(6.6) 



6.f 



M\uh,^ot)<mf\\L,^oi 



Now differentiate with respect to a;' G M. both sides of (6.6) m times to get 

(1)™^)* + (-f )™D"(^"'^i:i''i:i™M) + \D'^,u = D^'J/ 

in O^. Note that ^'"u satisfies (3.5). Thus by Theorem 3.1 again we have 



E P"i?>IL,(o+) < N\\f\\ 



i2(0+)- 



The lemma is proved. 



□ 



HIGHER ORDER SYSTEMS 



25 



Lemma 6.5. Let T G (—00,00] and A > 0. There exists N = N(d,n,m,S) such 
that, for u G W2' "(Oy) satisfying (3.5), 



provided that 



(6.8) 



in 0+. 



Proof. By multiplying both sides of the equation ( |6.8D from the left by D\"^u we 
get 

(6J) 

Note that 

5R(-l)'"(i:>2"u,wt)^+ = 1 /" |i:)™M|2(r,a;)da; > 0. (6.10) 
Indeed, this holds true because 



/ \D1'u\^{T,x)dx = / ^\D'^u\^dxdt = (i:>™u, i:>™ut)c,+ + {DTu,DY' 



and 



the latter of which follows from the boundary condition (5.5) and integration by 
parts. Hence by taking the real parts of (|6.9D and using (3.10) we have 



where a = (m, 0, • • • ,0). Thanks to the ellipticity condition and Young's inequality, 

2 

i2(0+) 



+Nie,6)Y,\\Dl"^-^Di, 



^(^)^'ll"llL(o+) + ^(^)ll/llL(oJ)- 



Choosing a sufficiently small e and using Lemma 6.4, we prove (3.7). 

Now we are ready to state and prove the main theorem of the section. 
Theorem 6.6. Let T G (—00, 00]. There exists N — N{d, n, to, S) such that 
E A1-^P"u|I^^(^+)<7V||u* + (-1)'"£ou + A«||^^(„J) 

for all X > and u G W2''^"^{0^) satisfying 



□ 



u{t, 0,x') = --- = D'^-^u{t, 0, a;') = 



on (—00, T) X 
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Proof. Thanks to Lemma 6.4 and interpolation inequalities, it suffices to prove that 



L2{0 



2m „ 



i2(0+) 



<iV|| 



Il2(oj;)' 



3.11) 



where f — Ut + (— l)™£oU + Xu. Lemma 3^ and Proposition |6.3| (with 2m in place 
of m) imply that 



This along with Lemma 6A and a sufficiently small e proves the inequality ( 3.11 ) 
without the ut term on the left-hand side. To complete the proof we simply note 
that 

ut = -{-irCou^\u + f. 

□ 



7. Mean oscillation estimates of some partial derivatives of 
solutions to systems on a half space 

The aim of this section is to derive several mean oscillation estimates of highest 
order derivatives of solutions to systems on a half space. Contrary to the whole 
space case, here we are only able to estimate parts of the highest order derivatives. 
More precisely, for divergence form systems, we give estimate of D^u, while for 
non-divergence form systems we present the estimate of D^u. We emphasize that 



these estimates alone are not sufficient for proving Theorem 21 and 2.6 
We still denote 



\a.\ — \j3\—m 



D°'A"^{t)D^u 



E 

|a I — — m 



A°'f>(t)D"'D^u. 



Recall that 



Q+{t,x) = Qrit,x)no+, Qt = Qrno+, 

Q; - (-r^", 0) X s;, = {x' e W'-^ ■. \x'\ < r}. 

7.1. Some auxiliary results for systems on a half space. We first prove some 
auxiliary estimates in this subsection. The first two are counterparts of Lemma 4.1 



and Corollary 4.2 



Lemma 7.1. Let <r < R < oo. Assume that u £ W^-'^'^ 



iQjl) satisfies 

u{t, 0, a;') = • • ■ = D'l'-^it, 0, x') = (7.1) 



on 



and 



in Q^, where f € L2{Q^). Then there exists a constant N = N{d,n,m,S) such 
that 



\D 



2m 



Furthermore, 

where N — N{d, n, m, S, r, R). 



Proof. By Theorem 6.£ the L2-estimate of s ystem s on a half spaces is available. 
Then the proof is the same as that of Lemma 4.1 with some minor changes. □ 
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Corollary 7.2. Let < r < R < oo. Assume that u G C'^^{Oto) satisfies (7.1) on 
Q'jzj and 



in Qjj. Then for any multi-indices 7 and d such that 

7 = (71,72, • • • ,7d), 71 < 2m, = (0,?92, • • • 

we have 

where N = N(d, n, m, S, r, _R, 7, i?). 



(7.2) 



Proof. From (7.2) it follows that 

D%t = -(-i)™^"'^i:'"i:'''i:'% 

in Q^. Each of the terms on the right-hand side is a constant times a term of the 
form D'^u, where I7I = 2m + \f3\ and 71 < 2m. Hence we only need to prove 

where 7 = (71, • • • ,7d) satisfies 71 < 2m. The proof of this inequality is identical 
to that of (4.7) in Corollary 4.2 with the only difference that, in I7I > 2m, we write 

where D'^u satisfies (7^) in as well as (7J) on Q'^, so that we can apply Lemma 

Next we derive a few Holder estimates of solutions. Throughout the paper, for 
a function g defined on a subset V in we set 

\git,x) - g{s,y)\ 



[g] 



sup ,K , I , 

{t,x),{s.^y)ev \t-s\2 + \x-y\'' 



where < < 1. 



Lemma 7.3. If u e C'^^{0^) satisfies (7J) on Q'^ and (7^) in , then 

MC1/2(Q+) <A^||u|li,(Q+). 

Proof. Let 

e+ = {(i,:Ei) G (-r2"\0) X (0,r)}, B,; = {x' £ M^^-^ : |a;'| < r}. 
The triangle inequality gives 



\u{t,x) - u{s,y)\ 



< 



sup 



\u{t,xi,x') - u{s,yi,x')\ 



\u{s,yi,x') -u{s,yi,y')\ 

sup — jjr^-jT^ ■.— ll+l2. 

^(s,yi)ee+ 



|a;'-y'|i/2 



To estimate Ii, we view u(t,xi,x') as a function of {t,xi) for a fixed x' e i3J. 
Then by the Sobolev embedding theorem 

\u{t,xi,x') - w(s,?/i,x')| 



<iV||H(.,a;')||^^i,.(e+) (7.3) 
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for each x' e B[. On the other hand, there exists a positive integer k such that, for 
each {t, xi) e Qi, 



sup \D{u{t,xi,x')\ + sup \ut{t,xi,x')\ 

j^Qx'£B[ x'eB[ 

2 

(7.4) 

where D-[u(t, xi, x') and Ut{t, Xi, x') are viewed as functions of x' G B'^. Combining 
(7.3) and (^) proves 



|7|<fc+2 
7l<2 



|i5|<fc 
i5i=0 



where the last inequahty is due to Corollary \!.2\ 

For the estimate of /2, we look at u{s,yi,x') as a function of x' S B[ for each 
{s,yi) G 8^. Again by the Sobolev embedding theorem, there exists a sufficiently 
large integer k such that 

\u{s,yi,x') - u{s,yi,y')\ 
,SVP , „ni/2 <N\Hs,yi,-)\\wiiB',)- 

x'¥=y' 

Moreover, as a function of {s,yi) G Qi, Dl.,u{s,yi,x'), j = 0, • • • , fc, satisfy 
fc fc 

for each x' £ B[. From the above two inequalities, we obtain (^^) with /2 in place 
of Ji . The lemma is proved. □ 

In the sequel, for a function g defined on O^, T G (— oo, cx)], we denote by £{g) 
{= £g) the even extension of g defined on Ot- 

Corollary 7.4. Let A > 0, = (to,0,a;o), where to G K and Xq G M!^^^ . Assume 
that u G Cf^^ioL) satisfies (^ on Q^(Xo) = (<o - 42™,to) x B'^{x'q) and 

ut + {-l)"" C^u + \u = {) (7.6) 
m (5^(Xo). T/ien t/iere exists N = N{d,n,m,6) such that 

m 
k=0 

Proof. By using a translation in t and a;', we assume that Xq — (0, 0). Let A = 0. 
In this case, the inequality in the corollary follows from 

[^>]cv.(q+)<^P'""IIl.(q+)- 

To prove this, we apply Lemma 7.3 to D"]u. This is indeed possible because D™u 
satisfies {7A_) on Q'^ an d (^.2 ) in . To prove the case A > 0, we follow the steps 
in the proof of Lemma 4.4. □ 
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Lemma 7.5. Let r G (0, oo), k G [64, oo), A > 0, and Xq ~ {to, xq) G Oto- Assume 
that u G C;'^^(0+ ) satisfies (O) on M x W^-'^ and m g+,(Xo). T/ien 



(If P!?^.) ~ (£(i?^^i))Q,^(X„)|)Q.(X„) + \H\£u - (5^/)q,(Xo)I)q.(Xo) 



(7.7) 



fc=0 



where N ^ N{d, n, m, 6) . 

Proof. Wc first prove that, using a scaling argument, it suffices to prove the in- 
equality (7.7) only for r = 16/k. Indeed, assume that the inequality (7.7) holds 
true for r = 16/k. For a given r G (0, oo), let rg — 16/k, R = r/rQ, and 



w{t,x) = u{R t,Rx). It is easy to check that w satisfies (|7.l|) on 



and 



wt + {-iy'A°'f^{R^'''t)D°'D'^w + XR^'^w = 



_ (7.8) 

in Q'^ro O^o), where Yq = (sq, ya) — {R^'^'^to, R^^xq) G O^. Then by the inequality 
(0) with r — ro applied to the system ( 7^ ) , we have 

\£iD^^w) - i£{D:]w))Q^^^y„ A + A^i?" (\£w - {£w)Q^JyJ 

^ QrQ{yo) ^ 

m 

< Nk-^-^\^-^^R'^-^{\£{D^w)\^)i 



Note that, for example, 

(f (i?'=H)Q.„(lo) - i?^f(Z?'«))Q.(X„). 

Thus the inequality ( |7.9| ) leads to the inequality ( [7.7| ) for arbitrary r G (0, oo). 
Now we assume r = 16/k. We consider two cases. 

Case 1: the first coordinate of xq > 1. In this case, we see that Q^r/ie^-^o) = 
Qur/iei^o) and u satisfies the assumptions in Lemma 4.5, especially, u satisfies ( [7.6[ ) 
in Q^r /lei-^o) ^nd u can be extended to a function in H^i^JM^^) without changing 
the values of u on Q^r /lei^o)- Hence by the inequality ( 4.12) wi th Qnr/iei^o) in 
place of QuriXo) (note that k/16 > 4), the left-hand side of (7/7) is controlled by 



3(^0)' 



which is less than the right-hand side of ( [7. 71 ). 

Case 2: the first coordinate of xq is in [0,1]. By denoting Yq = (to,0,Xg), we 
have 

Qr{Xo) C Q2{Yo) C Qs{Yo) C Q^r{Xo). 

By Corollary 7.4 applied to u with 2 and 8 in place of 1 and 4, respectively (this 
case can be seen using a scaling argument as above), we have 

{\£{Dl^u) - i£{D^.^u))Q^(^Xo)\)QAXo) < iVr^[f (i?;'}«)]ci/.(Q,(y„)) 
<iV^c-^^A^-^(|f(I?^)|^)|^(,.^)<iV.^-^5:A^-^(|£(i^M|^)|^^(^^). 



k=Q 



k=0 



The second term on the left-hand side of (7/7) are estimated similarly. 



□ 
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7.2. Mean oscillation estimates of for divergence type systems on a 
half space. Now we state and prove the main result of this section. 

Proposition 7.6. Let r e (0,oo), k £ [128, oo), A > 0, and Xq = (^0,2^0) G Ci- 
Assume that u G W^^iodOto) satisfies on M x W^-'^ and 

ut + {-irC^u + Xu= ^"/« (7-10) 

I Q I < rn 

in Q~^^{Xq), where fa € L2,ioc{0^), \a\ < m. Then we have 



mo^u) - i£{D^»)QAx„)\)QAx„)+>^H\£u-i£u)Q^^Xo)\)QAXo) 



k=0 



|a| <m 



(7.11) 



where N ~ N{d, n, m, 6) . 

Proof. Muhiplying u by an infinitely smooth function as ^ below, we see that 
( 7.10| ) can be extended to a system defined on 0+ without changing the values of 
u and fa on, for example, QK.r/2- Thus without loss of generality we assume that 
u e H^{0^), fa e ^2(0^), and (|]l^) is satisfied in 0+ . We consider only A > 0. 
Take a C G C^iQnr{Xo)) such that 

C - 1 on Q^rMXo), C = outside (to - (Kr)2",to + (Kr)^™) x B^rixo)- 



Let £q'^-' — A"^.D°'D^, where A"'? are the standard mollifications with respect to 



t of A°'^{t). Also let fa'' be infinitely differentiable functions approaching in 
^2(0^) as e \ 0. By Theorem 6T, there exists a unique solution w(^) e n^iO^), 



satisfying (|7.lD on 



, to the equation 



|q I <m 



in O^. Since fa^ and ^"j^ are infinitely differentiable, by the classical theory for 



higher order parabolic systems, w^*^^ is infinitely differentiable. Moreover, for any e, 



Thus by Lemma 7.5 (note that k/2 > 64) 

m 

<iV«-^^A^-^(|£(W^))r)|^^(^^). 



fc=0 



Set w^"') = u-'y(^). Then w(^) e H^iO+o) and it satisfies (7.1) on : 



and 



+ (-l)"4^)zi;(^) + Au;^^) = i?"(C/^^' + /a - ./l"^) + (-1)"(4'^ - A)^ 



in 0+ . Denote the right-hand side of the above equality by D"ga^ . We apply 
Theorem |6.1| to the above equation as one defined on Of^ so that we have 



i|<r 



1q|<t 
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In particular, 



L2(Q?(Xo)) 



<^ E ^^'''\\f'^^h.iQUXo))+I^'^' (7-12) 



I a I < rn 



k=0 



l..(Q^.(X„)) < ^ E ^^"'ll/l^^ll..(Q^.(Xo))+^^^^ (7-13) 

|a| <m 



for all sufficiently small where 

— 1/3|— m 

Note that, for the even extension £g of a function g defined on 0+ , we have 



|a|<m 



whenever Xq G Oi,. This combined with (7.12) and (7.13) gives 



\a\ <m 



g^'-*(|^<''*"'"')|')L™ 



<N A=-i(|£/W|')*,^,.v,, + M-'-*/<", 

Q|<m 



same inequahty with fa"^ in place of fa plus the error term 

(r-^-i + (Kr)-™-^)/(^). 
To finish the proof we let e \ 0. 



Now by following the corresponding steps in the proof of Theorem [4.6| we see that 
the left-hand side of the inequality ( 7.11 ) is less than the right-hand side of the 



□ 



Remark 7.7. Later we need to have the mean oscillation estimate (7.11) for all 

Xq G Ooo, instead of Xq £ Oi, for functions £{D'^u), £u, and £fa defined on Ooo 



if the equation ( 7.10 ) is satisfied in 0+ . In order to do this, in case Xq e OooXOtj, 
we let Yq be the reflection point of Xq with respect to the hyper -plane {{t,0,x') : 
f G M, a;' e R''"^}. By Proposition 7.6 we get the esti mate (7.11 ) with Yq in place 
of Xq. Then it is not difficult to see that the estimate ( 7.11 ) holds true as well for 
Xq using the evenness of functions involved. The same claim can be repeated for 
Corollary 7.8, Proposition 8.5, and Proposition B.6. 
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7.3. Mean oscillation estimates of D^Tu for non-divergence type systems 
on a half space. As a consequence of Proposition Z.f , we easily get 



Corollary 7.8. Let r e (0,oo), k G [64, oo), A > 0, and Xq = {to.xo) G Oto- 
Assume that u e W^;l"\0^) satisfies (0) onRx R'^-'^ and 

Ut + {-irCou + Xu = f 

in (5+,(Xo), where f g L2.ioc{0^). Then we have 



2 m 
fc=0 



where N ~ N(d, n, m, S) . 



Proof. Since satisfies (7.1) on 



Corollary 4.7 



^, we can proceed as in the proof of 

□ 



8. Estimates for systems with special coefficients on a half space 

The estimates in the previous section imply the iyp-estimate of D^u in the 
divergence case and that of D^V^u in the non-divergence case. In order to estimate 
the remaining highest order derivatives, by the interpolation inequality (Proposition 
S.3), it suffices to estimate D™u in the divergence case and Df^u in the non- 
divergence case. To this end, in this section we consider 



£ow = A{t)D( 



2m „ 



J=2 



D 



2m „ 



where A{t) = A""{t), d = (m, 0, • • • , 0). 

For this special operator, we have the following improved L2-estimate. 

Lemma 8.1. Assume that u e Cf^^{Oto) satisfies 

u{t, 0, x') = • • • = Dl'-^uit, 0, a;') = 

on Q'j^ and 

ut + (-l)™-Cou = 
in Q'^. Then, for any multi-index 7, we have 

P''^IIl.(q;^) + P""*IIl.(q+) < ^II"IIl.(q+). 
where N ^ N(d, n, m, 6, r, R, 7) . 



(8.1) 
(8.2) 
(8.3) 



Proof. As noted in the proof of Corollary |4.4 it suffices to estimate the first term 
on the left-hand side of (3^). Also, we only need to treat the case when the multi 
index 7 satisfies 7' = 0, where 7 = (71,7')- In fact, if the inequality ( ^.3| ) is shown 
to be true with 7' = and a smaller R, since D'^ u satisfies (8J) on Q'^ and ( ^.2[ ) in 
Q~^, we can replace u by D'' u in ( ^.3|) . Then the right-hand side, N\\D^ "IIl2(Q+)' 



is bounded by that of ( p.3[ ) by Corollary 7^2. Furthermore, by the interpolation 
inequality with respect to xi, it suffices to show 

l2/n 



(8.4) 
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for Z = 0, 1, 2, .... To prove the above inequality, we first observe that, thanks to 
(B^), we have 

d 



m+1 



Ut 



J=2 



in Q^. This together with ( ^.l| ) implies that (first with I = 0, then inductively) 



on Qfl- Moreover, Dl"^'u satisfies (S.2) in Q'^. Therefore, by Corollary 7.2 applied 
to Dl"^'-u we have 



2(i + l)m 



where r < ro < R. This implies (8.4) by an induction on I. 
As a consequence of the previous lemma, we get 



□ 



Lemma 8.2. Let u G Cf^c{Oto) satisfy (|^ on Q\ and (|J) in Qj. Then, for 
any multi-index 7, 



sup + sup \D^ut\ < N\\u\\ 



where N = N(d, n, m, S, 7). 



Qt 



L2(Qt 



Proof. This is deduced from Lemma B.l in the same way as Lemma 4.3 is deduced 
from Corollary |4.2|. □ 



Note that in the following Holder estimates the first inequality is for all D'^u, 
I7I — m, whereas the second inequality is for Di"^u only. Similarly we see D™u 
and D\"^u in the following lemma and propositions as well. 

Corollary 8.3. Let A > 0, Xq = {to,0,XQ), where to G M and Xq E 
that u e C^ciOt) satisfies (|t|) on Q'^{Xo) and 



I'' ^. Assume 



in Q^{Xo). Then there exists N — N{d,n,m,S) such that 



[f(i^™u)]ci(Q,(Xo))<A^E^"" 

fe=0 
2m 

fc=0 



\\£{D'^u)\U,iQ,iXo)), 



(Qi(^o)) 



|f(D'=u) 



\L2{Qi{Xo))- 



(8.5) 

(8.6) 
(8.7) 



Proof. Similar to the proof of Corollary [7.4| , we prove only the case A = and 
Xq = (0,0). As noted in th e pr oof of Lemma 3T, D\"^u satisfies (B.l) on Q'^ and 
(B.2) in . In this case, (3.7) follows immediately from Lemma B.2 applied to 

Lemma 



2 also shows that 

where the second inequality is due to the fact that u satisfies (8.1) and the boundary 
version of the Poincare inequality. This gives the inequality [iA). □ 
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Lemma 8.4. Let r G (0, oo), k G [64, oo), A > 0, and Xq = {Iq, xq) G Oto- Assume 
that u G C^^(0^) satisfies (O) on M x W^-'^ and on (3+ (Xq). Then 



fc=0 



k=0 



where N — N(d, n, m, S) . 



Proof. Thanks to the Holder estimates in Corohary 8.3, we process as in the proof 
of Lemma 7.1. □ 



From the above lemma, by following the steps in the proof of Proposition |7.6| we 
prove the following two propositions. 

Proposition 8.5. Let r G (0,oo), k G [128, oo), A > 0, and Xq — {to,xo) G Oj,. 
Assume that u G C^^^ci^^) satisfies (O) onRx R''-i and 



ut + (-ir£o^ + Au= D'^f^ 

|a| <m 

in Q~^^{Xo), where fa G L2joc{0^), \a\ < m. Then we have 

'm 
fc=0 



|Q;|<m 



QKr(Xo)' 



where N = N{d, n, m, 6) . 

Proposition 8.6. Let r G (0,oo), k G [128, go), A > 0, and Xq = {to,xo) G 



Assume that u G C;'^^(C'i) satisfies ( |8l] ) on M x M^^^^ anrf 

nt + (-l)"£ou + Au = / 
in (5+^(Xo), where f G L2,ioc(Ci)- ^/len we /lawe 

2m 
fc=0 



where N = N{d, n, m, S) . 
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9. Lp-ESTIMATES FOR SYSTEMS ON A HALF SPACE 

With the preparation in the previous two sections, we complete the proofs of 
Theorem 2.5 and Theorem 2.6 in this section. 



Proof of Theorem 2.5. RecaU that the leading coefficients satisfy Assumption 2.1 



(p). As before, we may assume that T — oo, p > 2, the lower-order coefhcients of 
jC are all zero, and u G C°°{Oto) vanishing on 0+ \ Qjig{Xi) for some Xi e Oi. 
In this case, it follows from Prop ositi on [z!^ (also see Remark [7.7|) and the proofs 
of Lemma 5.1 as well as Theorem 2.2 that 

in 



|a|<r 



k=Q 



|q I <m 



fe=0 

(9.1) 

for any ki > 128. 

Now we move all the spatial derivatives except DY'^{A°'°'D™u) to the right-hand 
side of dU), and add (-1)" Ei=2 ^i"" to both sides. Here d = (m, 0, • • ■ ,0). 
Then for any Qf^^riXo), K2 <E [128, oo), r e (0, oo), G and y e M+, we have 



ut + (-l)'"(I?™(A""(t,2/)7^™7.) +;Ei?2m^) ^ J2 D"f^ + {-iy'Y.^TD>^ 

3=2 \a\<rn J=2 

where fa = fa for \a\ < m, 

h^U- E {-ir\A''^{t,x)+Af^''{t,y))D^u+{-iy^A''"{t,y)-A^^{t,x))D^u, 



|/3|=m 



and 



/a -/a - E (-i)"^"^^^w + (-i)'"(^""(i,y)-^""(i,^))^' 



for \a\ = m, a ^ a. In the last two expressions, we used the fact that 



As a consequence of Proposition 8.5 and the proof of Lemma 5.1, for any K2 > 128, 
(|£:(D'"«) - (5(D'"z.))q^(;,„)|)q^(;,„) < 7V«.2-'E^^"^(l^(^'")l')i..(x„) 



fe=0 



|a| <T7i 



|a|— m,a7^^Q 
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Choose a /i £ (l,p/2). This estimate combined with the Fefferman-Stein theorem 
and the Hardy-Littlewood maximal function theorem gives 



m+i 



fc=0 



(9.2) 



From ( |9.2[ ) and Proposition f3J^, we get 
III?"?/ 



\a.\—m,a^a. 



fe=0 



Nk. 



|a| <m 



N{e)K'^^'\\DT,u 



(9.3) 



Combining (9.1) and ( |9.3| ) we obtain the desired estimate by first taking K2 suffi- 
ciently large, then s sufficiently small, ki sufficiently large, and finally p sufficiently 
small. □ 



Proof of Theorem 2.6. It suffices to establish the apriori estimate when T = oo, 
the lower-order coefficients of L are all zero, and u G C°"{Oto) vanishes on 0+ \ 
Qiig{Xi) for some Xi € Oto. We use the strategy in the proof of Theorem 2.5 and 

consider two cases. 

Case 1: p £ (2, oo). It follows from Corollary ^8 that 

) + ) < NnT^i ll/IL,(oi) 



LpiOto) 



(9.4) 



fe=0 



for any ki > 64. We move all the spatial derivatives except A'^°'D\"^u to the right- 
hand side o f (pT^ ), and add (—1)™ Yl'j=2 21,™" both sides. As a consequence of 
Proposition |8.6| and the proof of Lemma 5.1, for any K2 > 128, 



2m 



mor'u) - (£:(i?ru))Q„(^,) |)q^(^„) < nk^^ e ^^-^ (i^ 



k=0 



+Nk 



m+i 



+Nk. 



m+i 



|a[— 2m,a7^2a 

This estimate combined with the Fefferman-Stein theorem and the Hardy-Littlewood 
maximal function theorem gives 

2m 



/c=0 



m+i 



l-/'llip(0+) +^'*2 



E 



|i:>"w| 



Lp(0+)- 



(9.5) 



|a| — 2m,Q:7^2a 
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From (3^) and Proposition 3^, we get 



2m 

fe=0 
m+4 II „ii 



Lip (Ooo ) 



m+i 



{p^+e)\\Dl^ 



Nk. 



(9.6) 



Combining (9.4) and ( |9.6| ) we obtain the desired estimate by first taking K2 suffi- 
ciently large, then e sufficiently small, ki sufficiently large, and finally p sufficiently 
small. 

Case 2: p G (1,2]. Thanks to Case 1 and Remark 2.7, we already have the 
solvability of 

ut + {-!)"' Cou + Xu ^ f 

on the half space for any q G (2, oo) and A > 0. The same duality argument in the 
proof of Theorem |2.3| yields the solvability of the same equation for any q e (1, 2). 
We can repeat the argument in Section |8| to deduce a version of Proposition 
with 2 norms replaced by q norms. Inspecting the proof of Case 1, to finish the 
proof it remains to have a proper version of Corollary 7.8 with 2 norms replaced 
by q norms. 

We claim that Lemma 7.3 is still true with L2 replaced by ig, q e (1, 00), i.e., if 
u £ C^^(d^) satisfies (O) on and (Q in Q+, then 



MC1/2(Q+) < A^||"|li,(Q+). 



This easily yields the desired version of Corollary 7.8 by following the lines in 
Section [7|. H owever, the claim does not follow directly from the proof of Lemma 7.3 
because ( |7.4D doesn't hold if the Wj^'^ norm on the right-hand side is replace by the 
W^'^ norm when q is close to 1. To get around this, we use a bootstrap argument. 
We first note that under the assumption of Lemma 7.3, for any 1 < r < i? < 4, it 
holds that 



(9.7) 



This can be shown in the same way as Lemma 7.1 and i.l based on the global 
W^''^"^ estimate on the half space. By the Sobolev imbedding theorem and (^ 
we have 



I«IIl,,(q+) <^II"IIl,(q+) 



for any qi > q satisfying 



1 1 
— > - 

91 q 



1 



We iterate this bootstrap process for a finite many steps on a sequence of shrinking 
half cylinders, and get 



where qi > 2(c?+ 1). Now by the Sobolev imbedding theorem again, we deduce 

Ik|lcl/2(Q+) < A^|lu||i^(Q + ), 

which is exactly our claim. The theorem is proved. □ 
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Remark 9.1. From the bootstrap argument above, we actually can get a finer 
boundary estimate as follows. If u e ^qfdci^'^)^'! ^ (IjOo) satisfies (7.1) on Q'^ 
and (7^) in , then for any and e £ (0, 1), 

where = N(d, to, n, q, e). 

10. Systems on a bounded domain 



We present the proofs of Theorem |2.10| and |2.11| in this section. We first treat 
the non-divergence systems. In this case, the proof is quite standard by using 
the technique of flattening the boundary and a partition of the unity. We give a 
sketched proof for the sake of completeness. 



Proof of Theorem 2.11. First, in a same way as Lemma 4.1 by using Theorem 2.3 
instead of Theorem 3.2, we obtain the following interior estimate for any < r < 
R < oo, Qr C Qr C fir and A > Aq 

\Wt\\LAQr)+ E ^'"^P""IU.(Q.O <^ll/IU.(QH)+^ll"llip(Q«)- (10-1) 

\a\<2m 



Similarly, Theorem 2.6 yields a boundary estimate: let < r < i? < oo, / G 
Lp{ Q^), and p be the constant taken from Theorem 2.6. Then under Assumption 
Ml {p), for any A > Ao and u e W^-^"^{Q+), we have 

II^*IImq+)+ E ^'"^II^""IImq^)<^II/IImq+) + ^II"IIl,(q+)> (10-2) 

|a|<2r?j. 

provided that u = Diu = ... = D^C^^^u = on Q'^ and 

ut + {-!)"' Lu + Xu ^ f inQ+. 

It is well-known that the ellipticity condition is preserved under a change 
of variables. Take to G (— cxd, T), a point xq G dfl and a number rg = ro{fl), so that 

nnBroixo) = {x e Broixn) : xi > <j){x')} 

in some coordinate system. We now locally flatten the boundary of dfl by defining 

yi^xi- <j){x') := <l>\x), yj = Xj := j > 2. 

Under the assumptions of the theorem, $ is a (72m- i,i diffeomorphism in a neigh- 
borhood of Xq . It is easily seen that the leading coefficients of the new operator in 
the y-coordinates also satisfy Assumption 2A wit h a p ossibly different p. Thus, we 
can choose a sufficiently small p such that from ( 10. 2|) , for Xq = {to, xq) and some 
ri = ri{n) < ro. 



l'"t|Up(nTnQ,,(Xo)) + E 

\a \ <2m 



Lp{flTnQ^^{X„)) 



<mfKinrnQ.„iXo))+Nj2\\D^uh^ 



(OTnQ,.o(Xo))- 



(10.3) 



3=0 



Finally, a partition of the unity together with ( |10.1 ) and (10.3) completes the proof 
for a sufficiently large Aq. □ 
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Now we turn to the divergence case. We need to introduce a special mollification, 
which was used, for instance, in |3^ . 

Proof of Theorem 2.1(\ Again we only give an outline of the proof. The interior 
estimate is similar to that of the non-divergence case. Theorem 2J2 implies that, 
for any 0<r<_R<oo, QrC Qr C 51t and A > Aq, 



I a I < m 



|a| <m 



We also have the boundary estimate by Theorem 2.5: Let < r < i? < oo, / S 
Lp{ Q^), and p be the constant taken from Theorem 2.5. Then under Assumption 
(p), for any A > Aq and u G we have 



E A- 

I Q I < 7n 



|a| <m 
ym—1 ^ 



(10.4) 



provided that u — Diu = ... = D'^ u = on and 

ut + {-irCu + \u = f inQ+. 

Take S (— oo, T), a point Xq £ dfl and a number tq G (0, By Assumption 
2.9, locally in some coordinate system, we have 

nnBr„{xo) ^{xe Broixo) : x^ > (j)ix')}, 

and the local Lipschitz norm of is less than pi. The goal is to locally flatten the 
boundary of dQ. However, (f> is not smooth in this case since it is only assumed 
to be Lipschitz continuous. To construct a smooth diffeomorphism, we define a 
function cj) on R"^ by 



<^(.t) 



V{y')<f>{x' - xiy')dy'. 



and 



Here rj £ C^{B[) has unit integral. It is easy to check that (^(0,a;') — cj) 
\D^(j){x)\ < N{xiY~^ pi. We now define 

yi=xi- 4>{x) := $i(x), y^ = Xj := $^(a;), j > 2. 

As before, the leading coefficients of the new operator in the y-coordinates also 
satisfy Assumption with a possibly different p. After some straightforward 



calculations using ( 10.21 ) and Hardy's inequality, we conclude, for Xq = {to,xo) and 
some ri = ri{fl) G (0,ro), 

|a|<m |a|<m 



Npi J2 A' ^"^P""IUp(aTnQ.„(Xo))- 

\a\<.7n 



(10.5) 



Using a partition of the unity together with (10.4) and (10.5), we complete the 
proof of the theorem upon choosing a sufficiently large Aq and small pi . □ 
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11. Remarks on the ellipticity conditions 

In this section we discuss some other ehipticity conditions appeared in the htera- 
ture, and show how our resuhs can be extended to systems under those conditions. 

The following strong ellipticity condition has been widely used before; see, for 
example, |^. 

Assumption 11.1. For all £ M''+^ and complex vectors — |a| = 

m,i — 1, n, 

5ft I E U^'^<fit,x)]>S\e, (11.1) 
\\a\ = \P\=m J 

where 5 > 0. 

The next condition is called uniform parabolicity in the sense of Petrovskii, which 
has been used, for example, in |3^, |2^, |l^, We define a matrix- valued function 
on R'^+i X (R'*\{0}): 

q| — |/3|— m 

Assumption 11.2. Let Xj{t,x,(,), j — l,...,n, be the eigenvalues of A{t,x,^). 
Then, 

^{\,{t,x,0)>S, j = l,2,...,n, (11.2) 
for aU (t, x) e R'*+i and C e M'^ \ {0}, where S > 0. 

We still assume that all the coefficients are bounded and measurable. Clearly, the 
Legendre-Hadamard ellipticity condition (2.1) is weaker than the strong ellipticity 
condition. However, it is stronger than the uniform parabolicity in the sense of 
Petrovskii. 

11.1. The strong ellipticity condition. Since it is stronger than our assumption, 
all the results in th is pa per hold true under this condition. Moreover, we can take 
Ao = in Theorem 2.10 for divergence form parabolic systems without lower-order 
terms. In this case the solution u satisfies 

E < ^ E (11-3) 

|a| <m |q I <m 



Indeed, by the method of continuity it suffices to prove the estimate (11.3). Due to 
(11.1) and the Poincare inequality, we easily get the unique solvability for p = 2 as 
well as 

E < iV E ll/olU.(o.)- (11.4) 

|a|<m |Q|<m 

In the case when p > 2, we add (Aq -|- 1)u to both sides of the first equation of (2.7). 
By Theorem |2.10| , it holds that 

E < iVi E Wfo^hAnr) + NiMlL^n^y (11-5) 

|a|<m \a\<m 

Take pi G (p, oo) such that 1 — d/p > —d/pi. By Holder's inequality. Young's 
inequality and the Poincare-Sobolev inequality, we get for any e > 0, 

\\u\\L^{fiT) < N{e)\\u\\L^^nT) + 4u\\L^^{fiT) < N{s)\\u\\L^^nT} + N2e\\Du\\L^(^nT)- 
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Choosing e = l/(2iVi7V2) and using ( pTil ) and jlTll ), we obtain jlT^ ) for p > 2. 
The remaining case p S (1,2) follows from the standard duality argument. 

11.2. The uniform parabolicity condition in the sense of Petrovskii. As we 

noted, this assumption is weaker than the Legendre-Hadamard condition. Under 
this assumption, for the solvability of parabolic systems, we need to impose a 
stronger regularity assumption on the leading coefficient, that is, they are VMO in 
both t and x. More precisely, set 

osct,4A^P,Qr{t,x)) = -f \A"^{s,y)--f A"P\dyds, 

J Qr(t,x) J Q,.(t,x) 

and 

A* — sup sup sup osCx {A"^^ ,Qj.{t,x)^ . 

(t,2;)eR''+i r<i?, |Q| = |/3|=m 

Assumption 11.3 (p). There is a constant i?o G (0, 1] such that < p. 

Next we show that the results in Section || about parabolic systems in the whole 
space (Theorem 2.2 and 2.3) still hold true under the assumptions above. As a 
consequence, we obtain interior estimates for both divergence and non-divergence 
type parabolic systems. We note that, for non-divergence type parabolic systems, 
the corresponding interior estimate was established in a recent interesting paper 
|33[ (see Theorem 2.4 there) by using a completely different approach. 

By inspecting t he p roofs of the main theorems, it is apparent that if the L2- 
estimate Theorem |3.l| is proved for parabolic systems with constant coefficients 
under the uniform parabolicity condition, then the remaining arguments can be 
carried out as before with obvious modifications. Indeed, we have 

Theorem 11.4. Let T G (— cxo, cxo] and 

jCou= J2 D"{A'^^D^u), 

\a\ — \P\—m 

where A"^ are constants satisfying the uniform parabolicity condition (11.2). Then 
there exists N — N(d, n, m, S) such that, for any A > 0, 

^ X'-^\\D'^uh,^o^^<N J2 A^||/aI|L.(o,), (11.6) 

|a|<m |a|<7n 

ifue HTiOr), fc e L2{Ot), \a\ < m, and 



M( + (-i)"/:o" + Am= (11-7) 

I Q I < m 

in Ot- Furthermore, for A > and fa G L2{Ot), \ol\ < m, there exists a unique 
u G -H^^Ot) satisfying ( |lL7| ). 

Theorem 11.4 is probably known before. For example, it can be derived from the 
results in 35 ; see also Theorem 10.4 in Chapter VII of |Q. Instead of appealing to 
those general results, here we present a direct proof of it. We need an elementary 
lemma, which is verified by a direction computation. 

Lemma 11.5. Let 5 > and U he an n x n upper triangular complex matrix 
satisfying 

\U\<5-^, SftAj>(5, j = l,2,...,n. 
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where Xj , j ~ l,...,n, are the eigenvalues of U . Then there exist real constants 
£,Si > 0, depending only on n and S, such that for any x G C" 

'^(x^BUx) > Si\x\^, 



where B = diag{e' 



n-l 



, £,1} and x^ is the conjugate transpose of x. 



Proof of Theorem \L1.4 - It suffices to pr ove ( 11.6 ) when u e C^{Ot) and A > 0. 
We take the Fourier transform of ( |ll.7 ) in x and get 

ut+Aim^^'u + Xu^ (11-8) 

\a\<rn 

Let A(^) — Q^UQ be the Schur decomposition of A, where Q — Q{(,) is an n x n 
unitary matrix a nd U = J7(^) is an upper triangul ar m atrix. Let B be the diagonal 
matrix in Lemma 11.5 . Multiplying both sides of ( 11.8 ) by BQu and integrating 
on Ot give us 

{BQu, Qut)or + {BQu, UQ\(f"'u)or + HBQu, Qu)ot 

- E iQ"BQU,{iO''fa)or- (11-9) 

|Q|<m 

As in the proof of Theorem 

^{BQu,Qut)oT >0. 

By the Plancherel equality, 

X^{BQu,Qu)or > N{e)X\\u\\l^^Ory 

To estimate the second term of the left-hand side of ( 11. 9| ), we use Lemma 11.5 
and the Plancherel equality to get 

^{BQu,UQ\^\^"'u)or > Si{n,\e"'u)or > N{n,m,6)\\D"'u\\l^^^^y 
The real part of the right-hand side of (11.9) is bounded from above by 

N J2 P""IU.(o.)ll/olk.(o.) 



I a I < T71 I a I < m 



If "2 



L-2(Ot) 



for all e > 0. To complete the proof of ( 11.6 ) it suffices to use the interpolation 
inequalities and choose an appropriate e. □ 

Remark 11.6. In contrast, under Petrovskii's parabolicity condition, the Dirichlet 
boundary value problem of parabolic systems is in general not well-posed when 
d > 2, as pointed out in §10 Chapter VII of [|8|. However, in the case d — 1, 
relying on a linear transformation one can extend Theorem 11.4 to systems on the 
half space with the homogeneous Dirichlet boundary condition; see, for instance, 
§10 Chapter VII of Thus, all the results in Section || about systems on a half 
space or a bounded domain remain true in this case. 

Acknowledgement 

The authors are very grateful to Nicolai V. Krylov and the referees for helpful 
comments on the first version of the paper. 



HIGHER ORDER SYSTEMS 



43 



References 

Agmon S.: Lectures on elliptic boundary value problems. Prepared for publication by B. Prank 
Jones, Jr. with the assistance of George W. Batten, Jr. Revised edition of the 1965 original. 
AMS Chelsea Publishing, Providence, RI, 2010 

Agmon S., Douglis A., Nirenderg L.: Estimates near the boundary for solutions of elliptic 
partial differential equations satisfying general boundary conditions, I , Comm. Pure Appl. 
Math., 12, 623-727 (1959); II , ibid., 17, 35-92 (1964). 
[3] AusCHER P., Qafsaoui M.; Observations on W^'^ estimates for divergence elliptic equations 

with VMO coefficients, Boll. Unione Mat. Ital. Sez. B AHic. Ric. Mat. 5 (2002), 487-509. 
[4] Bramanti M., Cerutti M.: Wp'^ solvability for the Cauchy-Dirichlet problem for parabolic 
equations with VMO coefficients. Comm. Partial Differential Equations 18 (1993), no. 9-10, 
1735-1763. 

5] Byun S., Wang L.; Elliptic equations with BMO coefficients in Reifenberg domains, Comm. 

Pure Appl. Math. 57 (2004), no. 10, 1283-1310. 
6] Byun S.: Ifessian estimates in Orlicz spaces for fourth-order parabolic systems in non-smooth 

domains, J, Differential Equations 246, no. 9, 3518-3534 (2009). 
7] Chiarenza F., Frasca M., Longo p.: Interior VK^.p estimates for nondivergence elliptic 
equations with discontinuous coefficients, Ricerche Mat. 40 (1991), 149-168. 

8] : ly^'^'-solvability of the Dirichlet problem for nondivergence elliptic equations with 

VMO coefficients, Trans. Amer. Math. Soc. 336, no. 2, 841-853 (1993). 
9] Chiarenza E., Franciosi M., Frasca M.: L^-estimates for linear elliptic systems with dis- 
continuous coefficients, Atti Accad. Naz. Lincei CI. Sci. Fis. Mat. Natur. Rend. Lincei (9) 
Mat. Appl. 5 (1994), no. 1, 27-32. 

Dl Fazio C: estimates for divergence form elliptic equations with discontinuous coeffi- 
cients. (Italian summary) Boll. Un. Mat. Ital. A (7) 10, no. 2, 409-420 (1996). 

Dong H.: Solvability of second-order equations with hierarchically partially BMO coefficients, 
submitted (2009). 

Dong H., Kim D.: Elliptic equations in divergence form with partially BMO coefficients, 
Arch. Ration. Mech. Anal, 196 no. 1 (2010), 25-70. 
: Parabolic and elliptic systems with VMO coefficients. Methods Appl. Anal, 16 



[10 

[11 

[12 

[13 

[14 

[15; 

[16 
[17] 

[18 

[19 

[20 

[21 

[22; 

[23; 
[24 
[25; 
[26 



(2009), no. 3, 365-388. 



coefficients, Calc. Var. Partial Differential Equations, to appear (2010). 

Ebenfeld S., L-^-regularity theory of linear strongly elliptic Dirichlet systems of order 2m 
with minimal regularity in the coefficients. Quart. Appl. Math., 60, (2002), No. 3, 547-576, 

Ejdel'man S.D.: Parabolic Systems, (1969) North-Holland, Amsterdam, London. 

Friedman A.: Partial Differential Equations of Parabolic Type., Prentice-Hall, Englewood 
Cliffs, N.J, 2008. 

: Partial differential equations. Corrected reprint of the original edition. Robert E. 

Krieger Publishing Co., Huntington, N.Y., 1976. 

Giaquinta M.: Introduction to regularity theory for nonlinear elliptic systems. Lectures in 
Mathematics ETH Ziirich. Birkhauser Verlag, Basel, 1993. 

Gilbarc; D., Hormander L.: Intermediate Schauder estimates, Arch. Rational Mech. Anal, 
74(4) (1990), 297-318 

Haller-Dintelmann R., Hec;k H., Hieber, M.: LP-L'J-estimates for parabolic systems in 
non-divergence form with VMO coefficients, J. London Math. Soc, (2) 74 (3), 717C736 (2006). 

Kim D., Krylov N. V.: Elliptic differential equations with coefficients measurable with 
respect to one variable and VMO with respect to the others, SIAM J. Math. Anal. 39, no. 2, 
489-506 (2007). 

: Parabolic equations with measurable coefficients, Potential Anal. 26, no. 4, 345-361 

(2007). 

Krylov N. V.: Parabolic and elliptic equations with VMO coefficients, Comm. Partial Dif- 
ferential Equations 32, no. 1-3, 453-475 (2007). 

: Parabolic equations with VMO coefficients in spa<;es with mixed norms, J. Punct. 

Anal 250, no. 2, 521-558 (2007). 

: Lectures on elliptic and parabolic equations in Sobolev spaces, American Mathemat- 
ical Society, 2008. 



44 



H. DONG AND D. KIM 



[27] : Second-order elliptic equations with variably partially VMO coefficients, J. Fund. 

Anal. 257, 1695-1712 (2009). 

[28] Ladyzenskaja O. A., Solonnikov V. A., Ural'ceva N. N.: Linear and quasilinear equa- 
tions of parabolic type. American Mathematical Society, Providence, RI, 1967. 

[29] Leonard: S., Kottas J., Stara. J.; Holder regularity of the solutions of some classes of 
elliptic systems in convex nonsmooth domains, Nonlinear Anal., 60, 925-944 (2005). 

[30] Maz'ya v., Mitrea M., Shaposhnikova T.: The Dirichlet problem in Lipschitz domains for 
higher order elliptic systems with rough coefficients, preprint. 

[31] Maugeri a., Palagachev D., Softova L.: Elliptic and Parabolic Equations with Discon- 
tinuous Coefficients, Math. Res., vol. 109, Wiley— VCH, Berlin, 2000. 

[32] Miyazaki Y.: Higher order elliptic operators of divergence form in or Lipschitz domains, 
J. Differential Equations 230, no. 1, 174-195 (2006). 

[33] Palagachev D., Softova L.: A priori estimates and precise regularity for parabolic systems 
with discontinuous data, Discrete Contin. Dyn. Syst. 13 (3), 721-742 (2005). 

[34] : Precise regularity of solutions to elliptic systems with discontinuous data, Ricerche 

Mat. 54 (2005), no. 2, 631-639 (2006) 

[35] Solonnikov V. A.: On boundary value problems for linear parabolic systems of differential 
equations of general form, (Russian), Trudy Mat. Inst. Steklov. 83, 3-163 (1965); English 
translation: Proceedings of the Steklov Institute of Mathematics. No. 83 (1965): Boundary 
value problems of mathematical physics. III. Edited by O. A. Ladyzenskaja. Translated from 
the Russian by A. Jablonskii, American Mathematical Society, Providence, R.L 1967 iv-|-184 
pp. 

(H. Dong) Division of Applied Mathematics, Brown University, 182 George Street, 
Providence, RI 02912, USA 

E-mail address: Hongjie_DongSbrown.edu 

(D. Kim) Department of Applied Mathematics, Kyung Hee University, 1, Seochun- 
DONG, Gihung-gu, Yongin-si, Gyeonggi-do 446-701 Korea 
E-mail address: doyoonkimSkhu . ac. kr 



